A note on the stability of slip channel flows 
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We consider the influence of slip boundary conditions on the modal and non-modal stability of 
pressure-driven channel flows. In accordance with previous results by Gersting (1974) (Phys. Fluids, 
17) but in contradiction with the recent investigation of Chu (2004) (C.R. Mecanique, 332), we show 
that slip increases significantly the value of the critical Reynolds number for linear instability. The 
non-modal stability analysis however reveals that the slip has a very weak influence on the maximum 
transient energy growth of perturbations at subcritical Reynolds numbers. Slip boundary conditions 
are therefore not likely to have a significant effect on the transition to turbulence in channel flows. 



The advances in micro-fabrication techniques using 
polymeric or silicon-based materials has allowed to gain 
significant understanding on the behavior of fluids at 
small scales HHQ. One topic of current interest con- 
cerns the validity of the no-slip boundary condition for 
Newtonian liquids near solid surfaces 0, IE IE 0> HI • A 
large number of recent experiments on small scales with 
flow driven by pressure gradients, drainage, shear, or 
electric field have reported an apparent breakdown of the 
no-slip condition, with slip lengths possibly as large as 
microns. The slip length, A, is defined as the ratio of the 
surface velocity to the surface shear rate; A = corre- 
sponds to a no-slip condition, and A = oo to a perfectly 
slipping surface. 

Since the transition to turbulence in wall-bounded 
flows occurs at large values of the Reynolds number, 
studies in shear-flow instabilities have usually been out- 
side the realm of microfluidics. However, a set of recent 
investigations of the linear modal stability of pressure- 
driven flows in two-dimensional channels [3, flOl Hlj has 
reported that slip boundary conditions decrease the crit- 
ical Reynolds number, from Re = 5772 (its classical no- 
slip value obtained for Poiseuille flow) to Re « 100, in 
strong disagreement with early calculations of Gersting 
|l2j . Such results would potentially have a major impact 
on both turbulence and microfluidic studies. 

The goal of this note is twofold. First, we resolve the 
disagreement between the above cited results. A care- 
ful anal ysis of the derivation of the equations used in 
Refs. 0,ll3)ll3 reveals that incorrect slip boundary con- 
ditions on the perturbations were used in the modal sta- 
bility analysis. The use of the correct appropriate bound- 
ary conditions on the perturbations reveals the strongly 
stabilizing effect of slip on the eigenvalues of the linear 
stability operator, confirming earlier results Recent 
advances in the domain of shear flow instabilities have 
however revealed the usual lack of relevance of modal 
stability analysis, contrasted to non-modal stability anal- 
ysis, in subcritical transition in channel flows (for a re- 
view see, e.g., Refs. |l3|, UM)- The second goal of this 
note is therefore to quantify the effect of slip on the non- 
modal stability of viscous channel flows. To this end, we 



compute the maximum transient energy growth |15| in 
the presence of slip at subcritical Reynolds numbers. We 
find that, for all the considered combinations of stream- 
wise and spanwise wavenumbers, the effect of slip on the 
maximum energy growth and on the associated optimal 
perturbations is weak. 

Problem setting. We consider the flow between two 
parallel plates located at y* = ±h of a fluid with 
shear viscosity /i driven by a constant pressure gradient 
dp* /dx* in the ir*-direction. If we non-dimcnsionalizc 
lengths by h, velocities by U rc f = h 2 (—dp*/dx*)/2[i, time 
by h/U rc f and pressure by pU 2 ol , the dimensionless in- 
compressible Navier-Stokes equations for the velocity and 
pressure fields, (u,p), read 



— + u • V 



-VP + -W 
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where we have defined the Reynolds number for this flow 
as Re = phU xe {/fj,. We assume in this paper that the 
flow satisfies slip boundary conditions on both surfaces, 
with slip lengths Ai and A2 at y = h and y — —h, respec- 
tively. If we define the Knudsen numbers Kni = Xi/h 
and Kn2 = A2//1, and denote by (it, v, w) the streamwise, 
wall normal and spanwise components of u, the boundary 
conditions for Eq. JIJ are v = 0at?/ = ±l and 

du dw 
u + Km— = w + Km— = 0, y = l, (2a) 

du dw 
u-Kn 2 — = w-Kn 2 — = 0, y = -1. 2b 
dy dy 

Linear stability. We are interested in the stability of 
the steady unidirectional base flow U = U(y)e x satisfy- 
ing Eqs. GJ and 10, 



U(y) 



2(Kni + Kn 2 + 2KniKn 2 ) 



2 + Km 
2 (Km - Kn 2 ) 
2 + Km + Kn 2 



Kn 2 



y 



y 2 - 



(3) 
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In the absence of slip, Km = Kn 2 = 0, and Eq. J2J 
duces to the standard Poiseuille solution U(y) = 1 — 
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FIG. 1: Neutral curve uii(a, /?, Re) = for the symmetric slip 
case (a) and asymmetric slip case (b), (3 = 0, and values Kn = 
(no-slip), 0.01, 0.02, and 0.03. 

In order to characterize the stability of Eq. J3J), we write 
the total velocity field as the sum of the base flow plus 
small perturbations, u = U + u', p = P + p', and lin- 
earize the Navier-Stokes equations around (U, P ). T his 
procedure is classic and we refer, e.g., to Refs. 0, 
for the details. The same standard procedure is ap- 
plied to the boundary conditions (Eq. |2J . These linear 
boundary conditions are satisfied by both the total flow 
{u = U + u',v = v',w = w'} and the base flow itself 
{U, 0, 0}. Consequently, a simple subtraction shows that 
the boundary conditions for the perturbations are also of 
the form of Eq. [21 These boundary conditions are the 
same as those used by Gersting in his stability analysis 
|l2j and differ from the incorrect boundary conditions 
used in Refs. 0, 0, |n| that implicitly assume u' = 
at y = ±1. Therefore, in Refs. 0, 0, |n|, slip bound- 
ary conditions are assumed for the basic flow but no-slip 
boundary conditions are used for the perturbations, lead- 
ing to incorrect results. 

Following a standard procedure (see, e.g., Ref. 01 f° r 
details), the linearized Navier-Stokes equations are re- 
cast in a set of two differential equations for the wall- 
normal velocity, v' , and the wall-normal vorticity, rf — 
du' /dz — dw' jdx. Exploiting the homogeneous nature 
of the streamwise and spanwise directions, perturbations 
are Fourier-transformed in the form 

u'(x,y,z,i) = u(a,y,/M)e l(a ^ z) , (5) 

and therefore, r]'(x,y,z,t) = ff(a,y, (3,t) e l ( ax+ P z \ with 
rj = i[3u — iaw. The standard evolution equation for 
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FIG. 2: Critical Reynolds number for linear stability Re(Kn) 
for the symmetric and asymmetric slip cases. 

(v, rf) , is finally obtained to be [l4| 




where the operators are defined as 

C = -iaUA + iaD 2 U + A(A/Rc), (7) 

C = -i/3DU, (8) 

S = -iaU + A/Re, (9) 

with A = D 2 — a 2 ~ f3 2 , where D denotes derivatives 
with respect to y The fourth-order system of equations, 
Eq. ©, requires boundary conditions for both v and rj. 
Using the continuity equation, i au + Dv + i (3w = 
together with the boundary conditions in Eq. it is 
straightforward to show that the boundary conditions for 
(v, ff) are 

v = Dv + Kn^v = 0, y = l, (10a) 
v = Dv-Kn 2 D 2 v = 0, y = -1, (10b) 
lj+Kn 1 D?j=0, y=l, (10c) 
rf-Kn 2 Drf= 0, y = -1, (lOd) 

and for simplicity, we restrict the analysis in this note to 
symmetric slip (Kni = Kn2 = Kn) and asymmetric slip 
cases (Kni = Kn, Kn2 = 0). We emphasize again that 
these boundary conditions are different from those used 

m Refs. nnaim, where instead Kni and Kn2 were set 

to zero in Eq. I|10l) . 

Numerical method. A Chebyshev collocation method 
is used to discretize the system, Eq. and standard 
methods (described in Ref. [14J and references therein) 
are then employed to compute eingenvalues, eigenmodes 
and transient energy growth. The standard implementa- 
tion of these methods is modified by changing the stan- 
dard homogeneous no-slip boundary conditions into the 
more general slip boundary conditions (Eq. I10[l . All the 
results presented below have been obtained with 97 collo- 
cation points. Convergence of the results had been ver- 
ified and the code has been thoroughly tested by com- 
paring both the modal and the non-modal results in the 
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FIG. 3: Map of the iso- values of the transient energy growth 
G max (a, P) for Re = 1500 in two cases: No-slip (solid line) and 
symmetric slip boundary conditions with Kn = 0.03 (dashed 
line). The values of G max are 10, 100, 200, 300 and 400 from 
the outer to the inner curve. 



case of no-slip [F^ , as well as with the modal symmetric 
slip results reported in Ref. |T^ . 

Influence of slip on modal stability. The modal sta- 
bility analysis assumes solutions in the form of normal 
modes, {v, rf}(a, y, 0, t) = {v,rj}(a,y, 0,oS) e _Mt , where 
the complex frequency, u, is the solution to an eigen- 
value problem, which is solved numerically. The flow is 
found to be linearly unstable if there exists at least one 
eigenvalue with positive imaginary part, u>i > 0. The 
Squire theorem [1(3] applies to this flow and the critical 
modes are two-dimensional (i.e., with — 0). The neu- 
tral curve u>i(a,/3 = 0, Re) = in the symmetric slip 
case is displayed in Fig. Boundary slip is found to 
significantly shift the neutral curve towards larger values 
of the Reynolds number, indicating a strongly stabiliz- 
ing influence of slip on linear stability. Results for the 
asymmetric slip case, displayed in Fig. ^d, are similar, 
although less pronounced. The dependence of the crit- 
ical Reynolds number for linear stability, Re c , with the 
Knudsen number, Kn, is diplayed in Fig. [2] and confirms 
the stabilizing effect of slip on shear-flow instabilities. 
Our results, which use the correct boundary conditions, 
Eq. jlUII . agree with the symmetric slip calculations of 
Ref. |12(, but, as expected, are in strong contradiction 
with the conclusions reported in Ref. EH UM ■ 

Non-modal stability analysis. In the absence of slip at 
the walls, the Poiseuille flow is known to undergo tran- 
sition to turbulence at Reynolds numbers well below the 
critical Reynolds number corresponding to the onset of 
linear modal instability. This strongly subcritical transi- 
tion scenario has been related to the strongly non-normal 
nature of the linearized operator (Eq. [SJ) , explaining the 
potential of the flow to sustain large transient energy 
growth, possibly triggering the transition to turbulence 
for values of the Reynolds number much smaller than Re c 
[ilj. . The standard modal stability analysis is there- 
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FIG. 4: Evolution in time of the initial condition leading 
to the largest finite-time energy amplification G max for the 
Fourier mode (pt,f3) = (0, 2), Re = 1500, and in the symmet- 
ric slip case with Kn =0 (no-slip), 0.01, 0.02 and 0.03. 

fore extended to the non- modal (or generalized 01) sta- 
bility analysis where, for instance, the maximum tran- 
sient energy growth is computed. Let us define, for a 
given Fourier mode, the instantaneous kinetic energy of 
the flow perturbations as 

E(t,a,/3,u ) = J Jfi(a,y,/3,t)| 2 dy, (11) 

which is a function of time and the initial condition, Uo = 
u(a, y, j3, 0). If we denote by G{t) the energy growth at 
time t, maximized over all non-zero initial conditions, 

nu o\ { E(t,a,0,u o )\ 

G(t,a,0) = max — — — - , 12) 

then the maximum transient energy growth possible over 
all times, G max (ev,/3), is defined as 

G ma x(a,/3) = max G(t, a, 3). (13) 

In Fig. |31 we report the iso- values of G ma , x (a, 0) com- 
puted for Re = 1500 for both the no-slip (solid line) 
and the symmetric slip case (dashed line). Although the 
maximum energy growth with slip is always larger than 
in the case of no-slip, the increase is small and therefore 
slip hardly affects transient energy growth. The maxi- 
mum energy growth is obtained for a = and 8 = 2 for 
both slip and no-slip boundary conditions. Fig. 0] dis- 
plays the time evolution of the optimal energy growth 
G(t,a = 0, j3 = 2), at Re = 1500 and in the symmet- 
ric slip case, for different values of the Knudsen number. 
The small increase of the optimal growth with Kn ap- 
pears in all cases. Furthermore, the time where the max- 
imum growth is attained is also slightly increased by the 
slip. As both the square root of maximum growth and 
the time at which it is attained depend linearly on the 
Reynolds numbers, these effects suggest that the effect of 
slip induces an increase of an effective Reynolds number, 
which is consistent with the observation that slip flows 
have, for the same forcing, a larger flow rate than in the 
case of no-slip. 
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FIG. 5: Optimal initial condition, v opt (t = 0), and optimal 
response, r} op t(t = W), leading to largest transient energy 
growth at Re = 1500 and for (a,/3) = (0,2), for both no-slip 
(solid line) and symmmetric slip with Kn = 0.03 (dashed line) 
boundary conditions. 



the v' and rj variables, this means that the optimal 
initial perturbations are of v-type, with r\ negligible, 
while, on the contrary, the most amplified response is 
of Tj-type, with v negligible. This is also the case with 
slip boundary conditions. In Fig. [5] we reproduce, for 
Re = 1500, the optimal initial condition, v opt (y,t = 0) 
(left), and the optimal response f} op t(y, t = i max ) (right), 
corresponding to the largest transient energy growth, 
G max (a = 0,/3 = 2). The shape of the optimal initial 
perturbation differ slightly from the no-slip case, while 
the optimal responses are nearly undistiguishable, except 
near the wall, where the effect of the slip boundary con- 
ditions is apparent. The lift-up mechanism, by which 
low amplitude vortices are converted into large ampli- 
tude streaks seems therefore to be only slightly sensitive 
to slip boundary conditions at the wall. Similar results 
are obtained for other values of (Kn, Re) and for asym- 
metric slip boundary conditions. 



In the case of no-slip channel flow it is known that the 
initial perturbations inducing the largest energy growth 
are streamise vortices, while the most amplified response 
consist in streamwise streaks. Translated in terms of 
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